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IRTM I 1-pm channel brightness temperatures to also decrease in 
regions where low 20-pm channel brightness temperatures are 
observed [7,8]. The maps also show new phenomena, the most 
striking of which is a clear tendency for the low-brightness tempera- 
ture regions to occur at Fixed geographic locations. During this 
season, the coldest low brightness temperatures appear to be con- 
centrated in distinct regions, with spatial scales ranging from 50 to 
300 km. There are approximately a dozen of these concentrations, 
with the largest centered near the location of the south residual polar 
cap. Other concentrations are located at Cavi Angusti, and close to 
the craters Main, South, Lau, and Dana. Broader, less-intense 
regions appear to be well correlated with the boundaries of the south 
polar layered deposits, and the Mountains of Mitchell. We have thus 
far detected no evidence for horizontal motion of any of these 
regions. 

The fact that the low brightness temperature regions do not appear 
to move and are correlated with the locations of surface features 
suggests that they are not artifacts of the IRTM instrument or its 
viewing geometry, but the result of processes occurring on the sur- 
face or in the lower atmosphere. Presently, we do not know whether 
other low brightness temperature regions that have been observed 
during the southern winter or during the northern fall and winter 
exhibit similar spatial and temporal behavior. We intend to better 
understand the cause(s) and implications of these phenomena through 
modeling and further analysis of the Viking and Mariner 9 datasets. 
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Introduction: The near-surface martian wind and temperature 
regimes display striking similarities to terrestrial desert counter- 
parts [1,2]. The diurnal radiative cycle Is responsible for establish- 
ment of a pronounced thermal circulation in which downslope 
(katabatic) flows prevail during the nighttime hours and weak 
upslope (anabatic) conditions prevail during the daytime. The low- 
level wind regime appears to play an important role in modifying the 
surface of the polar regions [3]. Viking imagery of the north polar 
cap shows evidence of eolian characteristics such as dunes, frost 
streaks, and wind-scour features. The direction of the prevailing 
wind can in cases be inferred from the orientation of surface features 
such as frost streaks and ice grooves. 

For the past several years a numerical modeling study has been 
in progress to examine the sensitivity of thermally induced surface 
winds on Mars to the patterns of solar insolation and longwave 
radiative. The model used is a comprehensive atmospheric me- 
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Fig. 1 . Time evolution of (a) wind speed and (b) wind direction at the 
lowest sigma level (20 m) after 24-hr integration of constant slope runs for 
terrain slopes of 0.0005, 0.00 1, 0.002, 0.004, and 0.008, corresponding to 
curves A, B, C, D, and E respectively. 

soscale equation system that has been employed previously for 
simulation of Antarctic katabatic winds [4]. The model equations 
are written in terrain-following sigma coordinates to allow for 
irregular terrain [5]; prognostic equations include the flux forms of 
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Ffe. 2. Diurnal variation of(a) ground temperature, (b) wind speed, and(c) wind direction at the lowest sigma level for the solar cycle simulations over the 0.008 
slope for solar declination angles of 24', 1 2'. 0°. -1 2', and -24', corresponding to curves A. B. C, D. and E respectively. 


the horizontal momentum equations, temperature, and continuity. A 
surface-energy budget equation is also incorporated in which the^ 
surface temperature is determined. Explicit parameterization of 
both solar insolation and longwave radiation is included. Turbulent 
transfer of heat and momentum in the martian atmosphere is as- 
sumed to follow the similarity expressions in the surface boundary 
layer on Earth [6-8]. The vertical grid consists of 15 levels (o - 
0.998, 0.99, 0.98, 0.97, 0.96, 0.94, 0.92, 0.90, 0.85, 0.775, 0.70, 


0.60, 0.50, 0.30, 0.10). The high resolution in the lower atmosphere 
is necessary to capture details of the boundary layer flows. The 
lowest level corresponds to a height of approximately 20 m above 
the ground, the second level to 100 m. 

Probably the greatest uncertainty in model simulations of ther- 
mally driven martian flows is the specification of the terrain. While 
detailed topographic information is available for certain regions, 
only a fairly broadscale representation is possible over the martian 
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Fig. 3. Intensity of solar radiation (W nr 2 ) (a) reaching the surface and (b) corresponding to surface temperatures from numerical si mulation of diumal cycle 
over martian north polar cap. 


north polar cap. For three-dimensional simulations, terrain heights 
for the martian north polar region have been obtained from the U.S. 
Geological Survey map, which have been digitized onto a 57 x 57 
grid with a spacing of 75 km. 

Model Results: Two-dimensional simulations of thermally 

induced flows . A series of numerical experiments have been con- 
ducted that focus on the relationship between katabatic wind inten- 
sity and terrain slope. The model runs arc valid for high-latitude 
(75°), nocturnal conditions similar to midwinter on the north polar 
cap in which no solar radiation reaches the ground. A horizontal grid 
consisting of 20 points with a grid spacing of 20 km was used. The 
results of five 24-hour uniform slope runs are presented here. In each 
experiment, the model atmosphere was started from rest to isolate 
the katabatic wind. An initial lapse rate of 3°C knr 1 was used with 
a surface temperature of approximately 220 K at the vertical refer- 
ence level. 

In all cases, the katabatic wind reached a quasisteady state 
within the first 1 2 martian hours. The resulting evolution of the wind 
speed and wind direction at the lowest sigma level (approximately 
20 m above the surface) for the five martian katabatic wind simu- 
lations after 24 hr are shown in Fig. I. Curves A-E correspond to 
terrain slopes of 0.0005, 0.001, 0.002, 0.004, and 0.008 respec- 
tively. The Intensity of the katabatic wind (Fig. I a) appears sensitive 
to the terrain slope. For steeper slopes, the flow remains directed 
more downslope (downslope direction is 180°), which implies the 
horizontal pressure gradient force remains considerably greater 
than the Coriolis acceleration. 

To test the sensitivity of the martian slope flows to solar forcing, 
numerical experiments have been conducted in which the full cycle 
of solar forcing is replicated over sloping terrain. Results for a 
constant slope of 0.008 at 75° latitude will be described. 


Five numerical simulations have been conducted covering the 
seasonal range of solar declination angles (24°, 12°, 0°, -12°, and 
-24°); the model equations arc integrated for three complete mar- 
tian days to allow the model to settle into a stable diumal oscillation. 
The results presented here are taken from day 2; only minor varia- 
tions were seen beyond the first diumal cycle. It is assumed that the 
polar cap is composed of “dirty” ice with an albedo of 0.50. All 
simulations start from a rest state. Thus the influence of large-scale 
pressure gradients in the free atmosphere is neglected. This implies 
that all atmospheric motions arise due to the longwave radiative 
cooling or solar healing of the sloping terrain. No solar insolation is 
allowed for the first 12 hr of model integration time to allow the 
drainage flows to become established before model sunrise. 

Figure 2 illustrates the diumal course of the surface temperature, 
wind speed, and wind direction over the 0.008 slope martian terrain 
for solar declinations of 24°, 1 2°, 0°, - 1 2°, and -24°, corresponding 
to curves A, B, C, D, and E, respectively. The Sun never sets during 
the midsummer period and never rises for the winter case. The 
ground temperature (Fig. 2a) undergoes diumal oscillations of 30 K 
in summer, the magnitude of the oscillation decreases with the 
approach of the autumn; the diumal ground temperature oscillation 
amounts to 15 K at the equinox. Maximum temperatures appear an 
hour or so after local noon. Wind speeds at the first sigma level 
(Fig. 2b) show marked diumal trends during summer and equinox 
periods. Maximum wind speeds occur in the early morning hours 
coinciding with a minimum in the solar insolation in midsummer or 
just before sunrise in other simulations. Note that the simulated 
midsummer katabatic wind maximum of approximately 10.5 ms -1 
(reached in the early morning hours) is 3 ms -1 less than seen for the 
other cases. This reflects the insolation from the midnight Sun, 
which retards development of the katabatic wind. Wind directions 




38 Results from the MSATT Program 


streamlines 


SICma level = 2 



(b) viNO SPCCO <«•/*> SlCriA lOCL * 2 



Fig. 4. Mode! simulations of(a) streamlines and(b) wind speeds of slope flows at the second sigma level for diurnal cycle experiment. 


at the first sigma level throughout the diurnal course for the five 
numerical simulations are shown in Fig. 2c. The downslope direc- 
tion is 180° for these simulations. The Coriolis force acts to deflect 
the katabatic wind some 30° to the right of the fall line of the terrain 
for the winter katabatic wind case. The wind directions show sur- 
prisingly little variation with time except for the midsummer dec- 
lination angle of 24° (curve A). Note that upslope flow is modeled 
during the early afternoon hours of the summer case. The effect of 
solar insolation appears to retard but not overcome the katabatic 
forcing in all but the summertime case. This emphasizes the robust 
nature of the martian katabatic circulation. 

Three-dimensional simulations of the near-surface flows over 
the martian north polar cap . Numerical simulations have also 
been conducted incorporating solar insolation over the irregular 
terrain of the northern martian latitudes. For these simulations, the 
model equations have been integrated for a period of three martian 
days to ensure that the wind and temperature regimes have had 
ample time to adjust to the cycle of solar healing. It was found that 
the second- and third-day results of the model integration were 
nearly identical. Here the third-day results of one such experiment 
will be presented. In this experiment, the solar declination is as- 
sumed to be 25°, representative of midsummer conditions. It is 
assumed that the atmosphere is dust-free. As before, the model is 
initialized about a slate of rest to isolate the thermally induced slope 
flows in the lower atmosphere. Values of albedo and thermal inertia 
were taken from Paige [9]. 

Figure 3 illustrates the solar insolation received at the surface 
and corresponding surface temperatures. Local noon at this time is 
along the 60° meridian in Fig. I. Note that significant spatial 
variations in the intensity of solar radiation reaching the surface are 
found near the polar cap in response to albedo differences. Surface 


temperatures range from 265 K along the southern edge of the 
domain at local noon to around 190 K in the nighttime section of the 
model domain. 

Anabatic motions become established in response to the solar 
heating of the sloping terrain. Figure 4 illustrates the streamlines 
and wind speeds at the second sigma level for the time of the local 
solar forcing shown in Fig. 3. Note that especially along the local 
noon meridian, the streamline patterns are dramatically altered by 
the solar insolation as compared to pure katabatic flow. Upslope 
wind conditions are simulated at middle to high latitudes from 
morning through midaftemoon. As seen in the previous simulation, 
the strongest winds are associated with the steepest terrain. Upslope 
flows reach approximately 7 ms -1 in places along the local noon 
meridian. The intensity of the upslope flows appears to be weaker 
than the nocturnal katabatic winds. This is in agreement with 
numerous terrestrial observations such as over the Antarctic ice 
sheet. The katabatic wind speeds simulated over the nighttime 
sector centered are only slightly weaker than pure katabatic winds 
(not shown), emphasizing the rapid response of the wind field to the 
radiative- induced temperature changes near the surface. 

Summary: Numerical simulations of the martian near-surface 

wind regime using a mesoscale atmospheric model have shown that 
the thermally induced near-surface winds are analogous to terres- 
trial circulations. In particular, katabatic wind displays a striking 
similarity to flows observed over Antarctica. Introduction of solar 
radiation strongly perturbs the slope flows; anabatic conditions 
develop in middle to high latitudes during the daytime hours due to 
the solar heating of the sloping terrain. There appears to be a rapid 
transition from the katabatic to the anabatic flow regimes, empha- 
sizing the primary importance of radiative exchanges at the surface 
in specifying the horizontal pressure gradient force. 
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Over the past several years a number of investigators have 
described geomorphic evidence for and paleoclimatic significance 
of large standing bodies of water or ice sheets within the nonhem 
lowland plains of Mars [e.g., 1-14). The details of the timing, 
emplacement mechanisms, and sizes of these bodies differ mark- 
edly from one group of investigators to another, however. For 
example, Jons [ 1 ,2] envisioned a “mud ocean” covering much of the 
northern plains, with sediment slurries derived from a variety of 
peripheral sources, including the fretted terrains and outflow chan- 
nels. Lucchitta et al. [4] pictured an ice-covered ocean, fed by large 
circum-Chryse ice streams, analogous to those in Antarctica. Parker 
et al. [ 10, II ) indicated two or more highstands of a sea or ocean that, 
most recently, would have been charged by catastrophic floods, but 
may have existed more or less permanently during Noachian and 
Hesperian time. Interestingly, the shorelines of JOns’ “mud ocean,” 
Lucchitta et al.’s ice-covered ocean, and Parker et al.’s most recent 
sea, or “interior plains” [II], coincide almost precisely around the 
northern lowlands, though the details of the mechanisms by which 
key boundary morphologies are thought to have been produced 
differ. Baker et al. [ 1 2] pictured a plainswide ocean emplaced by the 
major outflow channels relatively late in martian history, and coined 
the term “Oceanus Borealis” for this ocean. 

Taking a more conservative approach to the question of standing 
water in the northern plains, Rotto and Tanaka [14] have relied on 
volume estimates of maximum discharge from the circum-Chryse 
outflow channels, which they feel limits any standing water to one 
or a few large, ephemeral lakes. The locations of these lakes are 
based on the identification of broad, shallow topographic basins on 
thepresent martian topographic maps [15]. Similarly, Scott et al. [ 1 3] 
have indicated evidence for several large lakes across the northern 
plains, some exhibiting connecting spillways, that were fed by a 
variety of channel sources peripheral to the plains. Delineation of 
these lakes is based on a similar assessment of the topography, but 
also included the identification of shore morphology. 

AH the above studies would have benefited greatly from the 
advent of the high-resolution topography afforded by the Mars 
Observer Laser Altimeter [16], which would have produced global 
topographic maps beginning in early 1994. For example, basin 
volume estimates in Parker et al. [11] are loosely based on the 
available topography with its very large vertical enors. These esti- 
mates, when compared to estimates by others of the water dis- 


charged by the Chryse outflow channels, suggest the possibility that 
the volumes required to fill the basin may be at or beyond the high 
end of the estimated volumes available from the channels. High- 
resolution topography is needed to sort out the common modifiers of 
shoreline elevation, such as tectonism, isostatic rebound, and sedi- 
ment desiccation and compaction, that probably altered the topog- 
raphy of the northern plains after the putative surface water was lost, 
so that the original topography can be reconstructed. Until a refiight 
of the laser altimeter or some similar instrument, elevations derived 
using the currently available high-resolution topographic tools — 
photoclinometry and shadow measurements — cannot be accurately 
tied to the global datum. It is not possible, therefore, to be certain 
that basin volumes based on the current global topography provide 
better than a crude approximation of the volume of ancient standing 
bodies of water in the northern plains. Until they do, such estimates 
cannot, by themselves, either point to nor preclude the presence of 
surface water or ice within the northern plains prior to the latest 
catastrophic floods; the uncertainties are still too large relative to the 
flood volume estimates. 

Can photoclinometry and shadow measurements be used to 
determine the volume of the basin without having to link the mea- 
surements to a global datum? Since the boundary, or shoreline of the 
basin cannot be tied to the datum and typically has no useful local 
relative height to measure, what is needed are a number of measure- 
ments of the height of the paleoshoreline(s) distributed across the 
basin — soundings, in effect. 

Parker et al. [9, 1 1 ] described a type of small knob in the northern 
plains that resembles terrestrial and lunar slcptoes (volcanic apron) 
and terrestrial wave-cut islands, and applied the nongenetic term 
stepped massifs” to the martian knobs. If these are upland outliers 
that had been abraded through wave action in an unfrozen ocean, or 
through ice-shoving in an ice-covered ocean, then the height of the 
apron above the surrounding plains could provide a measure of the 
local basin depth in the vicinity of the knob. Since stepped massifs 
are distributed over broad expanses in several places in the northern 
plains, it should be possible to measure the variation in basin depth 
regionally. With the exception of those regions where the available 
image scale is insufficiently small, it should be possible to measure 
the heights of the aprons to within a few tens of meters. 

As a feasibility test of this approach, photoclinomettic profiles 
are being compiled from Viking Orbiter images of the Cydonia 
Mensae region, which includes images with high Sun elevations 
(necessary to avoid shadows) and images with low Sun elevations 
(to enable the use of shadow measurements as an independent 
check) at high resolution ( 40-100 m/pixel). Both asymmetric and 
symmetric photoclinometic profile models are being used, and the 
results cross checked with one another to minimize the errors. An 
apron-height map, potentially a paleobathymetric map of part of the 
margin of “Oceanus Borealis,” will be compiled of this data to 
determine whether variations in apron height are consistent with a 
lacustrine interpretation. 
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